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Abstract 

The Boussinesq model of convection in a flat layer with heating from below is 
considered. We analyze the effects of anisotropy caused by rapid rotation in 
physical and wave spaces and demonstrate the suppression of energy transfer 
by rotation. We also examine the structure of the wave triangle in nonlinear 
interaction. The range of parameters is adapted to the models of convection 
in the geodynamo. 
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1. Introduction 

Approximation in the form of homogeneous and isotropic turbulence is 
quite crude for describing many geophysical applications. Usually, there are 
two reasons why this approximation can be violated. One of them is a rapid 
daily rotation characterized by a small Rossby number, Ro <C 1. This case 
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is typical for proble ms of meteo rology and physics of the ocean where its 



influence is principal ( Gil 



19821 ) . However, the most critical regime with an 



extremely small Rossby number occurs in the liquid planetary cores. There, 
the Coriolis force is three orders of magnitude larger than the nonlinear term, 
which changes the balance of the forces in the Navier-Stokes equation and the 
spectral properties of turbulence for large Reynolds numbers, Re. Another 
source of anisotropy of convection in the liquid core is a strong magnetic 
field which, in combi nation with rotation, leads to the plate-like structure of 



convective patterns (IBraginskv and Meytlisl . 



Donald and Roberts 



1990 



Matsushima et al. 



1999; 



20041 ). In the present paper we consider the former 



source of anisotropy related to rapid rotation. 

The introduction of rotation leads to substantial rearranging of the flow, 
both in the physical and in the wave spaces. In spite of this, the Coriolis force 
itself does not produce work. It can redistribute the energy between the scales 
and lead to inverse cascades, which is known in direct numerical simulations 
(DNS ) as an increase of the kinetic energy o n large sea 
predicted in the renormalization group theory (jMcCombl . 



cs f 


Hossain 


1992 


). As a 



19941 ) 



the slope of the spectrum of t he kinetic energy changes from —5/3 to —2 



((Zhou, 



1995 : 



Constantin 



2002). This ch ange is 



close ly related to the break 



of the energy transfer over the spectrum (jZhoul . Il995l ) . 

Let us proceed to the thermal convection problem, when an additional 
equation of heat transfer leads to the transformation of cellular convection to 



the cyclonic convection. This problem has been studiec 



oneering wo rks on thermal convection (1 Chandr asekhar 



thoroughly in t 



1961 



Busse 



19701 ). as well as in later papers on the geodynamo, see ( IRiidiger and Hollerbachl . 



Roberts 



re pi 



1965 



2 



2004 ) for references. 

It appears that rapid rotation causes a new kin d of balance b etween 



pressure and the Coriolis force (geostrophic balance) (jPedloskyl . 119871 ): l z x 

dV 

V ~ Vp, which leads to small gradients along the axis of rotation z: — — ~ 

n , 

q As a result, cyclones and anticyclones prolonged along z of a small 
diameter appear. If one takes the parameters of the Earth's liquid core and 
decreases the amplitude of the heating sources to the onset of convection, the 
diameter of cyclones will then be d c ~ E 1 ^ 3 L z fa 10~ 5 L Z , where E ~ 10~ 15 
is the Ekman number and L z is the height of the cyclone. In reality, the 
heat sources are quite larger, and there is a wave packet of such cyclones. 
However, estimates of the order of the energy indicates that, at least for the 
first 3-4 orders of the wave number k, geostrophic balance takes place. This 
exceeds the extent of geomagnetic spectra and should be taken into account 



in geodynamo models. This was already started in the deve 



tensor appro ach to anisotropic viscosi ty (IPhillips and Ivers 



opmen t of the full 



simulations (IMatsushima et al 



mass transfer (IDonald and Roberts . 



20031 ). numerical 



1999) and estimates of the anisqtropy of heat 



2004 



Matsui and Buffett 



2005|). 



The other role of the Coriolis force, which will be the subject of our paper, 
is to control the energy transfer over the spectrum in the system. Whereas 
if before the nonlinear term was of the same order with the pressure, which 
could not block the curl part of the term, then with the Coriolis force the 
pressure and other potential forces (e.g., the Archemedean force) block the 
potential part of the nonlinear term. And the rest of the curl parts of the 
forces and Coriolis force block the nonlinear term. As a result, the spectrum 



1 Note that balance of the potential part of the Coriolis force and the pressure holds. 
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becomes steeper. That is what we have on small scales, I <C d c . At the same 
time, on scales I 3> d c , the system can be in a state of statistical equilibrium: 
dissipation is negligible, and over the long term the energy exchange between 
the scales is small. It appears that the estimate of only the spectra of kinetic 
energy is not enough to resolve the states with and without rotation and 
intensive heat sources, and one needs to keep track of the fluxes of kinetic 
energy over the spectrum, which is considered below. 

Hereinafter, on the example of an 3D model of thermal convection in a 
plane layer, we consider the properties of anisotropy cause d by rotation both 



in the physical and especially in the wave spaces started in (IReshetnyak and Hejda . 



20081 ) typical for geodynamo regimes and compare these results with nonro- 



tating convection. 

2. The Boussinesq model 

2.1. Equations in physical space 

The thermal convection equations for an incompressible fluid (V • V = 0) 
in a layer of height L rotating with angular velocity Vt in a Cartesian system 
of coordinates (x, y, z) in its traditional dimensionless form can be expressed 
as follows: 



EPr" 1 



<9V 

^ + (V.V)V 



-VP - l z x V+ 

RaT^l z + EAV (!) 
^ + (V-V) (T + T ) = AT. 

Velocity V, pressure P and the typical diffusion time t are measured in units 
of k/L, pK 2 /L 2 and L 2 /k, respectively, where k is the thermal diffusivity, 



K V 

p is the density, Pr = — is the Prandtl number, E = — „ 



is the Ekman 



number, and v is the kinematic viscosity. Ra 



ag 6~TL . 



2Q.K 



is the modified 



Rayleigh number, a is the coefficient of vo l ume e xpansion, 5T is the unit 



20001 ) , g is the gravitational 



of temperature, for more details see (IJonesl . 

acceleration, and Tq = 1 — z is the heating from below. The problem is closed 

with periodical boundary con ditions in the fx, y) p lane. In the z-direction, 

j~ I Qy x 

we use simplified conditions ( Cattaneo et all . 120031 ): T = 0, V z = — - = 

oz 



dz 



at z = 0, 1. 



Equations in wave space 



To solve problem ([T]) we apply the pseudo-spectral approach (jOrszagi . 



197ll ). We follow the modification of the original approach that has been 
frequently used in thermal convection and geodynamo: we change FFT in 
z-direction to the sine-/cosine-decomposition leaving the full FFT in the 



horizontal directions. This let us t o provide the mentioned a 



conditions in z-direction , see also, (IJones and Roberts 



2000 



jove boundar y 



Buffett 



20031 ) . 



The equations are solved in the wave space. To calculate the non-linear terms 
one needs to apply the inverse Fourier transform, then calculate the product 
in physical space, apply the Fourier transform of the product, and finally 
calculate the derivatives in wave space. After eliminating the pressure using 
the divergence-free condition k • V = we arrive at: 



E 



Pr -i + k i Y 



dT 
~dt 



dt 



+ k 2 T 



J k 



(2) 



[(V- V)T + K] 



with 



• Pk : ~ k fc f k > k2 = k P k P> P = 1 • • • 3 

(3) 



Pr" 1 V x (V x V) + RaTl z - l z x V 

. k 

For integration in time we use the explicit Adams-Bashforth (AB2) scheme 
for non-linear terms. The linear terms are treated using the Crank-Nicolson 
(CN) scheme. To resolve the diffusion terms, we use the well-known trick 
of changing the variable which helps to increase the time step significantly. 
Consider equation 

w + eA = u - (4) 

Alter it to read 

8 -^l = U^ (5) 
at 

and then apply the CN scheme. 

The most time-consuming part of our MPI code are the Fast Fourier 
transforms. To make our code more efficient we use various modifications of 
known FFT algorithms, which take into account special kinds of symmetry 
of the fields. The optimal number of processors for the 128 3 grids is n ~ 50. 
Scalability tests have demonstrated even the presence of superacceleration if 
the number of processors < n. 

3. Fluxes in fc-space 

Now we introduce some diagnostic tools which are very helpful for the 
furthe r analy sis. To analyze the energy transfer in the wave space, we follow 



Frischl ( 119951 ). Let us decompose the physical field / into a sum of low- 
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frequency and high-frequency counterparts: /(r) = / < (r) + / > (r), where 



/<(r)= f* elkr > = E ^ e 

|fe|<x |fc|>if 



tkr 



(6) 



For any periodical functions / and g one has the relation (IFrischl . Il995l ): 



where 



df 9g 
ox ox 



< 9%- >= - < & >, <f > 9 < >=0, 
ox ox 



</(r) >=V~ 1 J f(r)dr 3 



(7) 



(8) 



stands for averaging / over volume V. Multiplying the Navier-Stokes equa- 
tion by V < one derives the equation for E < in sphere k < K: 



2- 1 EPr~ 1 



d V<V< 



Of 



+ U{K) 



RafT^V^) - E< (Vw < 



(9) 



where the integral flux of kinetic energy from the region k > K to k < K 
reads 

U(K) = (v<-(v r V 3 )v^, (10) 

with a summation over the repeating indexes i = 1 ... 3. Note, that flux of 
the Coriolis force is zero: -(V x V<) z = -(V< x V<) z - (V> x V<) z = 0. 
The first term in sum is zero due to vector identity and the second due to 
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(5). Introducing the local flux 



dU(k) 



T K (k) = / T K {k) dk = 0, (11) 



k=0 



where we have changed K to k, leads to the obvious relation for E in fc-space: 

dE{k) 



dt 

1 d 



T K (k) + F(k) + D(k), (12) 



where E{k) = < ^1 > ^ s the change of the kinetic energy at k, F(k) = 
Ra Pr d - 

— < TVf~ > is the work produced by the Archemedean force, and 

E dk 

D(k) = — Pi k 2 E(k) is the viscous dissipation. Equation (|12p describes the 
current flux of the kinetic energy through wave number k. 

To describe a triad interaction mechanism, we can also ask the question: 
"What is the form of the energy balance equation which describes the transfer 
of the energy from wave numbers Q and P to KT' A similar manipulation 
with the Navier-Stokes equation yields: 

^1 = T 3 + A(K) + D(K), (13) 
at 

where E(K), D(K) have the same form as in (pj), and T 3 = (Vt(K) ■ (Vj(P) ■ 
v/j^Q)}, A(K) = 5fp < T(K)V Z (K) >. 

It is also useful to introduce function T2(K,Q) = J T%(K,Q, P) dP. In 
general, one needs an explanation. It is possible to show that is the energy 



flux 



rom harmonic Q to harmonic K (see the references in ( lAlexakis et al. 



20071 )). The analysis of T 2 allows us to estimate whether the flux is lo- 
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cal or not, but not the locality of the interaction itself. The study of T 3 
allow us to recover the full structure of the wave triangle and to explain 
whether the interaction itself is local. We also note some useful properties of 
function T 2 : for the arbitrary periodical (or random homogeneous) nondiver - 



gent fields u(Q), w(fT)and V one has (jVerma 



2004 



Alexakis et al. 



20051 ): 



Tr(Q, K) = -Tr(K, Q), where T-HQ, K) = (u t {K) ■ (y 3 ■ V,)m t (Q) 
T WU (Q, K) = (wi(Q) ■ \ Vj ■ S/jjUi(K)j, which corresponds to the balance 
of the energy received by shell K from shell Q to that one given by Q to K. 
In the next paragraph, we consider the properties of fluxes Tk, T2, T3 on an 
example of model and find how they change when rotation is switched 



on. 



4. Basic properties of the fields 

The onset of convection in the plane infinite layer is a threshold phe- 
nomenon which occurs with the increase of the Rayleigh number to its crit- 
ical value, Ra cr . Here we consider three regimes of convection. We con sider 



simulations without rotation similar to (jMeneguzzi and Pouquet 



19891 ) (but 



without the magnetic field) and with rotation for two regimes with differ- 
ent amplitudes of the heat sources. All simulations were carried out using a 
quite rough grid, N = 64 3 , which helped us to present the statistics for the 
T3 fluxes. Here are the regimes: 

NR: Regime without rotation (the Coriolis term is dropped) , Ra = 9 • 10 5 , 



2 As there is no rotation, Ra and E do not retain their physical meaning defined in 
Section 2. 1. More details on parametriza tion of non-rotational magnetoconvection can be 



found in ( Meneguzzi and Pouquetl . 1989). 
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Pr = 1, E = 1, Re ~ 700. 
Rl: Regime with rotation, Ra = 4 ■ 10 2 , Pr = 1, E = 2 • 10~ 5 , Re ~ 200. 
R2: Regime with rotation, Ra = 1 • 10 3 , Pr = 1, E = 2 ■ 10~ 5 , Re ~ 10 4 . 

Regime NR in Fig{TJ corresponds to turbulent convection without rotation 
with a quasi-periodic in time behavior of the kinetic energy Ex{t). Note 
that a small-scale hydrodynamic helicity appears H n = V • rotV. For the 
case without rotation, because of the absence of a preferred direction, the 



me an helici t y Tiy is zero, see Fig|3l One can find the details of the problem 



in flMoffatt 



1978 



Krause et al. 




1980; 


Zeldovich et al. 


(Meneeuzzi and Pouquet. 


1989 


)■ 



19831 ). See also the 



Convection with rotation is characterized by the appearance of numerous 
vertical rotating columns (cyclones and anti cyclones). Thei r number is de- 



fined by the 



1965 



Busse 



Ikman number as k c ~ E 1//3 (I Chandr asekhar 



Jones 



1961 



Roberts! 



(12000) for details. For 



19701 ). see also the review by 
E 1 the number of cyclones and anticyclones is the same, whereas in the 
solar-like regimes, where the Coriolis and inert i al ter ms are comparable, the 



cyclones dominate (ISreenivasan and Davidson 



2008). For the liquid core of 



the Earth, E ~ 10 -15 . It leads to k c ~ 10 5 which is still imposs ible fo r DNS . 
Usually one is able to arrive at regimes with E = 10 -4 4- 10~ 6 ( Uonesl . 120001 ). 
The purpose of DNS is to find an asymptotic regime and then to extrapolate 
the results of the simulations to the terrestrial parameters. The side effect of 
such cyclonic convection is the decrease of the energy scale of the system and, 



3 As at the onset of convection the horizontal scales L x = L y of the convective cell are 
larger than the vertical scale L z , the box for simulations is usually longer in x-,y-directions. 
Here L x = L v = 5L Z = 5L. 
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as a result, the increase of viscous dissipation: the critical Rayleigh number 
depends on the Ekman number as Ra cr ~ E -1 ^ 3 . 

Regime Rl corresponds to the geostrophic state near the onset of convec- 
tion, see Fig{TJ The increase of Ra (regime R2) leads to the suppression of the 
regular cyclones and appearance of the small scale flows in the z-direction, 
and deviation from the geostrophic state to the so-called quasi-geostrophic 
turbulence, see Figj2j The nonlinear term tends to the amplitude of the Cori- 
olis force and pressure gradient, the temporal behavior of the system becomes 
chaotic. In both the cases (Rl, R2) the mean non-zero helicity generates: 
U H =< V ■ rotV > xy (the average is taken over the (x, y) plane), see Figj3j 
As was mentioned before, TC is zero for regime NR. For R2, H n is close to 
the linear function in the main volume. 

It is instructive to analyze the behavior of the kinetic energy as a function 
of depth for different components of the velocity field, see FigJH The different 
behavior of the transverse E% = (V^ + V y 2 ) /2 and longitudinal e\ = 
components of the energy depends on the different boundary conditions. 



For the nonrotating case NR, the values of the components are comparable: 

Ei? 

~ 1, where the factor 2 appears because of the summation 

z=0.5 



^NR — iT 

'II 



2E K 

over the two horizontal directions. In contrast, for regime Rl in the vicinity 
of the central part of the volume, the transverse velocity changes its sign 
and, as a result, the ratio increases: Ari ~ 0.1. This reflects the symmetry 
property of the flow at the threshold of generation. It is interesting that 
the increase of Ra (regime R2) leads to reinforcing the anisotropy in the 
whole volume in the opposite direction: Ari ~ 5, i.e. at the rather large 
Ra and under strong rotation, vertical motions are suppressed with rotation. 
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This is the so-called degen eration of t 



iree dimensional convection to the two- 



dimensional predicted by (IBatchelor 



19531). see also the s imilar degeneration 



Kraichnan 



19651 ). Also, the increase 



caused by a strong mean magnetic field 
of Re (R2) goes with the appearance of the layer structures at z — 0, 1. The 
latter is connected with the formation of a thermal boundary layer, thickness 
8 T ~ Ra~ 1/3 . 

5. Spectra of the fields and fluxes in the wave space 



Here we continue our analysis of the fluxes in k-space performed in lReshetnyak and Hejda 



( 120081 ) and consider anisotropy effects. For our regimes, we shell consider 



the integral spectrum of kinetic energy Ex{k), as well as its longitudinal 



4(*ii) = / / 



E K (k x , k y , k z ) dk x dk y and transverse E K (k±) = J E K (k x , k y , k z ) dk z 
spectra, where k\ = k 2 x + ky, k\\ = k z . For regime NR, the spectral estimates 
are close to the Kolmogorov dependence ~ k~ 5 ^ 3 , see Figj5l In other words, 
the spectra for NR are isotropic. It corresponds to the isotropic form of the 
convective cells at the main scale of the box. 

The spectra of convection with rotation differ from those for NR. Regime 
Rl is close to the onset of convection. The integral spectrum demonstrates a 
well-pronounced maximum, which corresponds to the scale of cyclones k c ~ 8. 
The increase of Ra (R2) fills the gap for k < k c and the spectrum starts to 
resemble the spectrum without rotation. 

Spectra E% and e\ for Rl are different: spectrum e\ does not feel 
boundary k = k c , whereas E^ at k < k c is close to the white noise, and 
at large k± the spectrum decays. For larger Ra, E\ tends to Kolmogorov's 
asymptotic, and the transverse spectrum is still white for small k and close 
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to the one without rotation at k± > k c . 

Note that for Rl and R2 the form of the integral spectrum Ek is defined 
by the form of its transverse component E K - 

The apparent similarity in the spectra for NR and R2 does not mean 
that the physics of the processes are similar. In this respect we remind the 
reader that the two-dimensional turbulence with inverse cascades and the 
three-dimensional turbulence with dire ct cascades have the same s l ope o f 
"-5/3" for the kinetic energy spectrum (IKraichnan and Montgomery! . Il980l ). 
However, the directions of the energy transfer through the spectrum in these 
two examples are opposite. 

Figj6] shows the fluxes of the kinetic Tk energies for the regimes mentioned 
above. At first we consider the integral over all directions. Regime NR for Tk 
demonstrates the well-known behavior for the direct Kolmogorov cascade in 
3D. For large scales Tk < 0, these scales are donors and provide energy to the 
system. On the other hand, the harmonics with large k absorb energy. The 
two-dimensional turbule nce exhibits mirror-symmetrical be havior relative to 



19801 ) . In this case the 



the axis of the abscissa (IKraichnan and Montgomery! , 
energy cascade is inverse. 

Rotation essentially changes the behavior of the fluxes of kinetic energy. 
The leading order wave number is k c . For k > k c we also observe the direct 
cascade of energy Tk > 0. The maximum of Tk is shifted relative to the 
maximum of the energy to large k; the larger Re, the stronger the shift. 
For k < k c , the behavior is more complex: for small k, the inverse cascade of 
kinetic energy takes place, Tk > 0. It is very tempting t o associate this small 



region with the appearance of coherent structures, see (iTabelingi . 



20021). On 
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the other hand, for the larger region of k (0 . . . k c ) we still have the direct 
cascade Tk < 0. The increase of Re leads to the narrowing of the region 
with the inverse cascade and to the increase of the inverse flux. One may 
suggest that the change of the sign of flux T K at k < k c is connected with 
the appearance of the non-local energy transfer: so that the energ y to the 
large -scales ki comes from modes |k 2 | ~ |k 3 | ^> |ki|, ki = k 2 + k 3 (IWaleffd . 



19921 ). Hence, in the case of rotation, two cascades of kinetic energy (direct 



and inverse) take place simultaneously. As show n in the simulations for the 



higher resolution (iReshetnyak and Hejdal . 120081 ) these results do not depend 
on the presence of the magnetic field, if the kinetic and magnetic energies 
are comparable in order of magnitude. 

To estimate the anisotropy of the fluxes, we decompose Tk into the sum 
of the transverse and longitudinal T K parts, so that in Eq. lflO]) the sum- 
mation over index j was used for j = 1 ... 2 and j = 3, respectively. For 
all three regimes, the fluxes in the z-direction are much smaller than in the 

horizontal plane, however, the form of T K are similar to that of Tk- This 

dE 

asymmetry can be explained as follows. Obviously, flux 11^ ~ — — is 
small for geostrophic flows. For NR, it is also small because the kinetic 
energy does not change too much along the length of the plumes in the z- 
direction in the main volume. As regards the boundaries z = 0, 1, V z is small 
there. That is why the main contribution to Uk yields the transverse flux: 
11^ ~ • V±Ek- This asymmetry distinguishes all three regimes from the 
homogeneous, isotropic models. 
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6. Locality of the energy transfer 

Let us consider the structure of the triad interactions. Figf7]is a diagram 

of the antisymmetric (with respect to the diagonal K = Q) fluxes T 2 for 

regimes NR, Rl, R2. For NR, the results are similar to the case with the 
I ^ Jil 

Alexakis et al.l (120071 ): harmonics with K > Q take the 



imposed force in 



energy from harmonics with K < Q (direct cascade of the energy). The 
maximum of the energy flux corresponds to the closest to the diagonal har- 
monics with K ~ Q, i.e. the local energy transfer exists. Note that there are 
regions (i.e. Q = 5, K = 20), in which there is a nonlocal inverse cascade of 
energy. Omitting these details, the whole behavior is close to the idealized 
Kolmogorov scenario. As a matter of convenience, we present the diagrams 
as a function of K — Q in FigJHJ This clearly demonstrates the existence 
of the direct cascade, as well as the local interactions with the local energy 
transfer. 

The rotation changes the behavior of T 2 at k < k c , leaving it unchanged 
at k > k c . We shall consider it in more detail. The positions of the extrema 
are near to k c , where the energy transfers to the region of k larger than the 
leading mode k c . On the other hand, at k < k c the flux decreases, which 
corresponds to the approach to the state of statistical equilibrium, observed 
in Figj6j On the larger scale (see FigJ7]^3)) we resolve the finer structure at 
small k: region A with a direct cascade of energy, but with an equiprobable 
energy transfer from the small Q ~ K, as well as from the rather large 
Q ~ 4K. Region B with an inverse cascade (as well as in region A with the 
weak energy flux of the order of 1/10 of its maximal absolute value) has an 
elongated strip-like form from Q ~ K to Q ~ This corresponds to the 
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occurrence of the small negative minimum at K > Q q 

The increase of Ra (regime R2) leads to the shift of the region with the 
inverse cascade in the direction of the small Q and K > Q. As before, we 
observe the equilibrium state at k < k c . This regime exhibits a longer interval 
of the wave numbers at k > k c with local energy transfer and the direct 
cascade. The relative input of the region with the inverse cascade becomes 
smaller, see FigJH](3) and at the same time this input shifts to the large-scale 
field at k -C k c , which can be interesting for geodynamo applications, where 
k c ~ 10 5 and the geomagnetic field generation regions (k ~ 1 -r- 10 3 for the 
typical magnetic Reynolds number R m ~ 10^ ~r~ 10 3 ) are separated at least 
by some orders of magnitude of k. 

7. Locality of interactions 

Let us consider the properties of function T3 for three regimes in more de- 
tails. Because of the symmetry of the problem, we expect that T 3 (K, P, Q) = 
T 3 (K, Q,P), which was used in constructing the discrete analog of the oper- 
ator. The case without rotation, see Figj9]^l, 2) demonstrates a quite inter- 
esting result: the largest input to the energy flux for mode K comes from 
two sources: from the P < K,Q <C P and Q < K, P <C Q- In other words 
the wave triangle (K,P,Q) is isosceles with a small angle between the equal 
sides Z(K,P) or Z(K,Q) . Taking into account that, according to FigJHl 
the energy transfers to K from the closest smaller wave number, we come to 
conclusion that the third small vector is a catalyzer in the interaction while 

4 We do not consider the behavior of function T2 at K < Q (regions C and D), because 
of its antisymmetry. 
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participating in the interaction this high-frequency wave does not provide to 
the wave K itself with energy. This result is beyond the classical Kolmogorov 
scenario, according which the wave triangle is an equilateral, i.e. not only the 
energy transfer is local, but the interaction as well. There are two reasons 
for this discrepancy. Having in mind the results of 
connect this result with a shortness of the considered spectra. On the other 
hand, it can be caused by anisotropy concerned with the gravity force. In 
this connection, it is known that, even for the larger Reynolds numbers in 



Alexakis et al 



(120071 ). we 



the problem with t he imposed externa 



can be substantial (jZhou. Yeungi . 



1996|). 



force, the anisotropy on small scales 



Our modeling demonstrates that there is a strong correlation of the fluxes 
on small scales with a buoyancy force on main scale. For this purpose, 
we introduce function r(K) = J T3 • max(P/Q, Q/P)dPdQ/ J T^dPdQ 
(P, Q > 0), see FigHTJ. Virtually, for all K the angle 

a(K) = J f-T 3 dPdQ/ J T 3 dPdQ (f = (P 2 +Q 2 -K 2 ) /2PQ) between vec- 
tors p and q is about 100°, i.e. the modes, which take part in the interaction, 
are orthogonal in the wave space. 

In the case with rotation the structure of the triangle with small Ra differs 
from the case without rotation. The particular wave interacts with the set 
of the waves - see the cross-structure of the diagram in FigJTOl This can be 
done by increasing the angle between vectors p and q. Note that P ~ Q for 
small K and P , see FigJTOl 

There is a well-pronounced inverse cascade of energy for small K, see 
FigfjJJ harmonic K = 7 receives energy from the higher wave numbers, trans- 
fering the energy, in its turn, to shorter waves. For the larger K, there is a 
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state with P ~ Q and a ~ 100°. In contrast to the case without rotation, 
the energy input to K is made by waves P ~ Q m 0.7K. The increase of Ra 
leads to the intermediate state. 

8. Conclusion 

To summarize the main results of the paper: the rotation sufficiently 
changes the structure of the flows in the physical space and its spectral 
properties. The introduction of rotation leads to the transformation of the 
typical Rayleigh-Benar cellular convection to cyclonic form. In general, the 
rotation suppresses convection due to enhanced dissipation: the horizontal 
scale of the cyclone is E~ 1//3 times smaller than its vertical scale L z . Moreover, 
the rotation leads to the selective suppression of the motion along the axis 
of rotation (Ar 2 <S 1) and it also upsets the mirror-reflecting symmetry of 
the system, which is reflected in the non-zero mean hydrodynamic helicity 

W± o. 

The behavior of the system in the regime with rotation with moderate 
Ra is very different for k > k c and k < k c . For k < k c there is a weak 
inverse cascade with nonlocal interactions. It is possible to speak of statistical 
equilibrium when the energy exchange between the Fourier modes is absent. 
For k > k c , the cascade is direct, however, the nonlocal flux from small k ~ k c 
exists. The slopes of the spectra for k < k c and k > k c are different. For long 
waves, the spectrum is close to white noise, and for large k, the spectrum 
decays as A; -3 . 

We have demonstrated that, even in a pure hydrodynamic system without 
a magnetic field, a variety of different interactions occurs in the Fourier space. 
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This analysis could be used for adjusting semi-empirical turbulent models 
proposed for geodynamo purposes in the future. In such models, there should 
be agreement of the spectra and fluxes in the wave space for scales larger than 
the cut-off scale d a with DNS of finer resolution. The nontrivial point here 
is the reproduction of the inverse cascades with d c < d a . 
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Figure 1: Sections of the vertical component of velocity, V z , without rotation (regime 
NR, two upper plots) at x = 4.3, and z = 0.8. Ranges of the filed: (—257, 506) and 
(—254, 572). Two lower plots demonstrate sections of V z with rotation (regime Rl) at 
x = 4.3 and z = 0.8. Ranges of the filed: (-88, 127) and (-55, 86). 




Figure 2: Sections of the vertical component of velocity V z with rotation (regime R2) at 
x = 4.3, and z = 0.8. Ranges of the filed: (-407, 449) and (-456, 600). 




Figure 3: The profiles of the hydrodynamic helicity H H for E = 3 • fO 5 , Pr = 1, 
Ra = 4 • 10 2 (solid line) and E = 3 • 1(T 5 , Pr = 1, Ra = 1.2 • 10 3 (squares) (1). Values are 
normalized at the mean over the volume kinetic energy (Ek) (2). Circles correspond to 
the regime NR. 




0.2 ■ 

J ■ ■ ■ ■ ■ 

0.2 0.4 0.6 0.8 1 




z 



Figure 4: Dependence of the full kinetic energy Ek (1), transverse E K (2) and longitudinal 
Ek (3) components for regimes NR (solid line), Rl (crosses) and R2 (diamonds). All 
quantities are normalized at Ek- 




Figure 5: Spectra of kinetic energy Ek, E^, as a function of k (1), k± (2) and fc| (3) 
components for regimes NR (solid line), Rl (crosses) and R2 (circles). The straight line 
corresponds to the Kolmogorov's law of ~ fc~ 5 / 3 . 




Figure 7: The fluxes of kinetic energy T 2 (k) for regimes NR (1), Rl (2 and 3), R2 (4). 



10' 



^ 



10 f 




■30 -20 -10 10 20 30 

K-Q 



10' 



^ 



10 



9 




■30 -20 -10 10 20 30 

K-Q 



4 10 11 • 






• 




V 


\ 




-4 10 11 • 







■30 -20 -10 10 20 30 
K-Q 



Figure 8: The fluxes of kinetic energy T 2 (K - Q) for regimes NR (1), Rl (2), R2 (3). 
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Figure 9: Fluxes T 3 for fixed K. NR: K = 5 (graph 1), K = 15 (2), Rl: K = 7 (3), 
K = 15 (4), R2: if = 5 (5), K = 15 (6) normalized at the maximal values for all P, Q 
and K. 




Figure 10: Dependence maximum of the ratio of two wave numbers P/Q, Q/P (left plot) 
and angle between vectors p and q (right plot), as a function of the resulting harmonic K 
for regimes NR (solid line), Rl (crosses) and R2 (diamonds). 



